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We implement the impedance measurement technique (IMT) for characterization of 
interferometer-type superconducting qubits. In the framework of this method, the interferome- 
ter loop is inductively coupled to a high-quality tank circuit. We show that the IMT is a powerful 
tool to study a response of externally controlled two-level system to different types of excitations. 
Conclusive information about qubits is obtained from the read-out of the tank properties. 



PACS numbers: 74.50.-|-r , 85.25.-j 
I. INTRODUCTION 



Quantum effects in mesoscopic superconducting cir- 
cuits of small Josephson junctions have attracted re- 
newed attention. It was clearly demonstrated that 
Josephson devices can behave like single microscopic par- 
ticles if they are sufficiently isolated from the environ- 
ment. Therefore, ideas developed in atomic and molecu- 
lar physics can be used for description of artificially fab- 
ricated circuits of macroscopic size. These concepts are 
stimulated further by the perspectives of a promising way 
to realize quantum bits (qubits) for quantum information 
processing. 

Qubits arc two-level quantum systems with externally 
controlled parameters. Generally, two kinds of such de- 
vices with small-size Josephson junctions have been de- 
veloped. One approach is based on charge degree of free- 
dom, basis states of this kind of qubits are distinguished 
by the number of Cooper-pairs on a specially designed 
island. The alternative realization utilizes the phase of 
a Josephson junction (or the flux in a ring geometry), 
which is conjugate to the charge degree of freedom. Due 
to macroscopic size of superconducting qubits, they are 
extremely sensitive to external disturbances. Thus, a 
backaction of a detector should be as small as possible. A 
lot of different detectors have been suggested in literature 
(see Ref. Q and references therein) . 

In this paper we review our results obtained on su- 
perconducting qubits by the impedance measurement 
technique (IMT). Below we shall discuss several quan- 
tum effects including macroscopic quantum tunneling, 
Landau-Zener transitions, Rabi oscillations, and direct 
resonant spectroscopy of the qubit energy levels. Finally, 
we present our very recent results of investigation of two 
coupled qubits. 



II. MACROSCOPIC QUANTUM TUNNELING 

For the flux qubits the Josephson energy dominates 
over the charge energy, Ej ^ Ec- It was predicted, that 
such systems should exhibit various quantum-mechanical 
effects including macroscopic quantum tunneling (MQT) 
of the flux.^ Indeed, predicted effects had been observed 
experimentallyi^i^i^ In this section we briefly discuss the 
main properties of the flux qubits and demonstrate that 
the IMT technique is a powerful tool for the investigation 
of the MQT. 

One of the realizations of the flux qubit is a super- 
conducting loop with low inductance Lq, including three 
Josephson junctions (a 3JJ qubit). ^ Its total Josephson 
energy is Ej = Ya=i Eji{(pi), where and Eji = HIci/2e 
are the phase difference and Josephson energy of the ith 
junction, respectively. Due to flux quantization, only 
(/)i.2 are independent with 03 = —(pi — 02 — Stt^o/'&o 
(<f>o — h/2e is the flux quantum) for negligibly small Lq 
(though sea^). 

At $0 = ^^oj the potential [/(0i,02) has two shal- 
low minima. These two minima correspond to the qubit 
states il^L and ipu, carrying opposite but equal supercur- 
rents around the loop. Therefore, according to the laws 
of quantum mechanics, near degeneracy the system can 
tunnel between the two potential minima. 

In the basis {tpL,tpR} and near $o = ^^o^ the flux 
qubit can be described by the Hamiltonian 

A £ 

H = -^f^x- ^(Jz- (1) 

At bias £ = 0, the two lowest levels of the qubit an- 
ticross (Fig. QJl) with energy gap A. With e changing 
sign, the qubit can either adiabatically transform from 
to V'fl staying in the ground state E- or switch to the 
excited state E+. The probability of the latter (called 
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FIG. 1: (a) Quantum energy levels of the 3JJ flux qubit vs 
external flux. The dashed lines correspond to the classical 
potential minima. In all graphs, the states A, B, C corre- 
spond to, say, '^L with left-rotating supercurrent. As $e is 
increased, these lose classical stability in favour of the corre- 
sponding states 'Pfl, denoted by D, E, and F. (b) Internal vs 
external qubit flux. 

Landau-Zener transition) for linear sweep e(t) — vt 
and £ changing from — oo to oo was calculatedSi to be 
Plz = exp(-7rA72?iz^). 

In order to demonstrate the principle of the IMT mea- 
surements of this system, let us consider the internal 
flux representation (Fig. ^) instead of the energy one 
(Fig. . A similar picture is usually used for the ex- 
planation of the operation of the conventional radio- 
frequency (rf) SQUID. The main difference between an 
rf SQUID and a qubit behavior is the existence of the 
adiabatic trajectory BE for latter one (see Fig. b). 
Let us assume that BE trajectory is forbidden and 
"qubit" is inductively coupled to the high-quality res- 
onant circuit. Then the system exhibits hysteretic be- 
havior >iS The tank circuit is simultaneously driven by a 
dc bias current Idc and an ac current I^f of a frequency 
Lo close to the resonance frequency of the tank circuit. 
Both currents produce the total magnetic flux applied 
to the qubit $e = '^dc + '^rfcosut. If the amphtude 
$r/ > ^ft) where is the half- width of the hystere- 
sis loop AC ED (Fig. the tank circuit will register 
the energy losses proportional to the loop long 
as |$(ic — ^*&ol < '^rf — ^h- These losses occur due 
to the jumps from E+ to E- at the ends of the loop. 
This idea was used by A. Silver and J. Zimmerman to 
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FIG. 2: Electron micrograph of the qubit at the center of the 
tank coil. 

build the first rf SQUID magnetometers.^^ If > 
the rf voltage across the tank circuit is a ^o-periodical 
function of applied dc flux Vri^dc) with local minima at 
= i<I>o -I- n$o, where n is integer. 
Now, let us take into account the additional "quan- 
tum" trajectory BE (see Fig. ^p). If its probability 
1 — Pi^z is nonzero but less than 1, two new closed 
paths ABED and BCEE are possible. There are two 
contributing trajectories, adiabatic and Landau-Zener 
transition. Therefore the net dissipation is Pioss = 
2-Plz(1 — Plz) and vanishes if Plz is either too small or 
too largei^. Due to the exponential dependence of Plz 
on the sweep rate, in practice this makes the quantum 
losses observable only if the bias sweep narrowly over- 
shoots the anticrossing, i.e., if 

l$<ic-i$o| (2) 

when <I>e changes slowly. Plotting Vri^dc) for $r/ > 
^h/'2,, a plateau flanked by two peaks is expected. The 
position of the dips depends on $r/ as follows from 
Eq. 13 Therefore in contrast to Vxi^dc) dependence of 
an rf SQUID the qubit should exhibit two local minima 
(in one period) which are symmetrical with respect to 
^dc = For the amplitude $r/ > the ACFD 

hysteresis becomes closed as well. Here, similar to the rf 
SQUID, on the K/(^'dc) dependence should appear the 
local minimum exactly at = \^o- Note, that <I>e here 
plays a role of bias e for the Hamiltonian . 

To test the ideas discussed above, we prepared litho- 
graphically square-shaped Nb pancake coils with induc- 
tance Lt on oxidized Si substrates for the tank circuits. 
We used an external capacitance Ct to be able to change 
the resonant frequency lot = \/^/LtCt- The line width 
of the 30 coil windings was 2 /im, with a 2 /im spac- 
ing. The quality factor of the tank was Qt ~ 1500 
at LOT ^ 20 MHz. The 3JJ qubit structure was fabri- 
cated out of Al in the middle of the coil by conventional 
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FIG. 3: Tank voltage vs. magnetic flux bias near the degen- 
eracy point of the qubit ^$o- From the lower to upper curve, 
the driving voltage is 10.2, 10.7, 11.2, 13.1, 17.2, 21.3, 22.0 
fiV rms (data vertically shifted for clarity). 



shadow evaporation technique. The Josephson junctions 
with critical current density jc ~ 300A/cm^ have areas 
« 130 X 620 nm2, 120 x 600 nm^, and 110 x 610 nm^ re- 
spectively. The loop area was 90 /im^, with Lq — 39 pH. 
The fabricated structure is shown in Fig. |21 

We measured Vri^dc) by a three-stage cryogenic am- 
plifier placed at « 2 K}^ Results for small driving voltage 
are shown in Fig. |3| For the smallest voltages no dissi- 
pative response is observed; the two "quantum" peaks 
appear around 10.7 /iVji^ and subsequently move apart 
without significant broadening. The "classical" peak ap- 
pears in the center, and with an ac bias threshold double 
the one of the quantum peaks — both as predicted above. 

Now assume that the probability of a Landau-Zener 
transition is small and the qubit adiabatically changes 
from ijjL to -0/?,, always staying in the ground state 
This means that the hysteresis ACFD is "shunted" by 
the BE trajectory. Therefore there are no losses caused 
by the flux jumps in the qubit. However, in the vicin- 
ity of B (see Fig. ^) the small change of the exter- 
nal magnetic flux causes a considerable change of the 
flux inside the qubit. Due to coupling of the qubit to 
the tank, the effective inductance of the tank-qubit sys- 
tem is changed, which leads to the change of the res- 
onant frequency. In this context the convenient mea- 
sure of that change is the imaginary part of the total 
impedancei^ expressed through the phase angle x be- 
tween driving current Ibiasit) = lac cos ujt, and tank volt- 
age Vrit) = Vt cos(wt -|- x)- For small Lg and if the am- 
plitude of /,/ is negligible, the results are summarized 



tanx = k^QxL, 



(3) 
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FIG. 4: Tank phase shift vs. flux bias near degeneracy 
fx ~ 0. Dots correspond to experimental data, solid line 
is the theoretical fit with A/h = 650 MHz. 



efficient, with M being the mutual inductance between 
them. The ground-state curvature isi^ 



$g(4i?2A2/J + A2)3/2 ' 



where 



f =^-- 
$0 2' 



(4) 



(5) 



and A is the normalized flux-to-energy conversion factor. 
Since all quantities in Eqs. ©"© can be measured in- 
dependently, experimental results can be compared with 
theoretical expectationSii& 

Strictly speaking, Eq. (|3Jl describes a measurement of 
the quantum object with vanishing back-action. There- 
fore, its validity should be provedJii Taking into account 
the influence of the tank on the qubit, the Hamiltonian 
(0) should be rewritten as: 



H 



e 
2 



-'J,{Qo + f + -fiT) + HT + HgB, (6) 



where 7 = IqM is the coupling coefficient between the 



qubit 's current, Iq 



and the current in the tank 



where k = M / ^JUqU^ is the tank-qubit coupling 



It- An infinitesimally small auxiliary force j{t) is re- 
quired for calculations of qubit 's magnetic susceptibil- 
ity. A heat bath operator Qq ^^^id a Hamiltonian ElqB 
describe internal mechanisms of dissipation and fluctua- 
tions in the superconducting loop. The high-quality tank 
is treated here as a quantum cavity is characterized by 
creation/annihilation photon operators a"*", a which obey 
the Bose commutation rules [a,a'^]_ = 1 etc. Quantum- 
mechanical operators of a current and a voltage in the 
tank are defined as: It — yhwTl^Xrr{a + 0+),!^ = 
iy^fiZ7r/2C^(a+ — a). For the Hamiltonian of the tank 
driven by a bias current luas and coupled to its own heat 
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bath Qf, we get the expression 

—LrlThias + Htb- (7) 

The internal heat bath of the tank Q;,, characterized by 
a free Haniiltonian Htb, results in a finite life time of 
the photons, 7^^, and, because of this, in a finite quality 
factor, Qt — lut/it- Assuming that h — l,kB = 1, we 
derive the Heisenberg equations for the tank operators: 

It = Vt/Lt, and 



1T-J^+<^T j^T ^ S.b + XuJt^z + -^ibias, (8) 



angle x- It follows from the averaged Eq. Hll() that am- 
plitudes of harmonic oscillations of the tank voltage and 
the bias current are related through 



Vre 



It 



LiT 
X2, ,2 

■ —j—XzzK^) 



^ (12) 



with Xzz(w) and x"z(w) being the real and imaginary 
parts of the qubit magnetic susceptibility ifTUIl . In the 
case when the tank is driven exactly with the resonant 
frequency, uj ~ lot, the voltage amplitude, Vt, can be 
found from the equation: 



where £,T{i) is a fluctuation source with zero average 
value, (^) = 0, and a correlator, {^b{i^)^b), which is pro- 
portional to the linewidth of the tank jt and the tank 
temperature T: {^b{t^)£,b) = {2-fTT / Ct)uj'^ ■ Because of in- 
ductive coupling the current and the voltage in the tank. 
It and Vt, affect the qubit current, Iq = IqUz- Using the 
linear response theory this influence can be presented as 
follows 



A 



(9) 



Vt 



'^{[eLqiiu^T^^A^T]? 

+ [jT + k'Lql'qUJTX'L{^T)?} 



-1/2 



(13) 



whereas for the voltage-current phase shift we obtain the 
expression 



tanx = -2k^LqI^QTx'zz{^T) 



(14) 



Here Qt = wt/(7t) is an effective quality factor of the 
tank wherein a broadening of tank's line width. 



where an operator az.o describes fluctuations of the qubit 
current caused by its internal heat bath, Qo, which is 
not correlated with the heat bath of the tank, Qb- We 
also take into account relations {S/SIt) — X{S/Sf) and 

It -~ Vt I Lt- 

The function {5<7;,{t)/Sf{ti)) involved in Eq. © is pro- 
portional to the derivative of the qubit current Iq{t) = 
{Iq) over the flux $y = LtIt created by the tank, 
6Iq{t) /S^{ti), or to the second derivative of the qubit 
energy profile, i?($), over the flux, d'^E{^)/d^^ (com- 
pare with Eq. O. It is convenient to characterize the 
qubit response on the action of the tank by means of the 
magnetic susceptibility-'^^ Xzz{^) defined as 



Then, the voltage in the tank obeys the equation 



(10) 



dti 



dt^ 



A2 



dt 



VTih) 



£_b + XioipOzS) 



1 

Ct 



X^'r-T^-Ibias- 



(11) 



It is evident from this equation that the tank voltage 
contains information about the magnetic susceptibility 
Xzziy) of the qubit. Similar to classical case this infor- 
mation can be extracted from measurements of the phase 



7T = 7T + k^Lqllu^TXzzi^T) 



(15) 



due to the qubit is taken into account. The magnetic 
susceptibility of the qubit (Eq. I10|l is calculated from the 
Bloch equations written in the form: 



X I 



(16) 



where and Ty are qubit's damping rates, ct" = 
— (A/wc) tanh(wc/2T) is a steady-state polarization of 
the qubit with energy splitting Wc — VA^ + e^, which 
is much higher than the resonant frequency of the tank, 
ujc ^ liJt- Because of this the decoherence and relaxation 
rates drop out of the expression for the magnetic suscep- 
tibility: 



Xzz (wt) = x'zzi^r) = 
A2 



'(A2 +e2)3/2 



tanh 



2T 



(17) 



As a result, the phase angle between the voltage in the 
tank and the bias current is given by the formula 



tanx 



A 



A2 



3/2 



X tanh 



VA2 



2T 



(18) 
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By making use of simple algebra it can be shown that 
at T = Eqs. ||2J) and IjlHIl are equivalent. Therefore, 
indeed measuring tan x as a function of the bias applied 
to the qubit, let us to determine the qubit's tunneling 
rate A. 

In order to realize the adiabatic response of the qubit 
experimentally, we fabricated 3JJ Al qubit with the fol- 
lowing parameters. The area of two, nominally equiv- 
alent junctions was about 190x650 nm^ while one is 
smaller, so that a = Ej3/Eji^2 ~ 0.8. The value of 
the critical current for larger junction was determined 
to be Ic ~ 380 nA. Qubit inductance, tank parameters, 
and measurement setup were the same as in the case of 
Landau-Zener transitions described above. 

The measured xifx) curve at nominal mixing-chamber 
temperature T = 10 mK is shown in Fig. ^ The curve 
was fitted by Eq. J^J) with A as a free parameter. Calcu- 
lated curve for the best fit parameter A/h = 650 MHz is 
also shown Fig. ^ This value of energy gap was in good 
agreement with the gap determined independently from 
temperature measurements^ 



III. RABI SPECTROSCOPY 

Quite generally, a two-level quantum system (includ- 
ing qubits), should exhibit coherent (Rabi) oscillations in 
time in the presence of resonant irradiation. Microwaves 
in resonance with the spacing between qubit's energy lev- 
els will cause their occupation probabilities to oscillate, 
with a frequency proportional to the microwave ampli- 
tude. Indeed, coherent oscillations between energy levels 
of the superconducting qubit were detected.--^'^*''^^'^^ 

In this section we show that the tank can be used for 
detection of Rabi oscillations as well. If the resonant 
microwave signal is applied, the phase-coherent oscilla- 
tions of the level occupation will only last for a finite 
time, which is usually called the coherence time. The 
correlation between the occupations can be expressed by 
an autocorrelation function or its Fourier transform, the 
spectral density. For example for the IMT, when the flux 
qubit is coupled inductively to a tank circuit, the spectral 
density of the tank-voltage fluctuations rises above the 
background noise when the qubit's Rabi frequency fipj 
coincides with the tank's resonant frequency wt- This 
forms the basis for our measurement technique of Rabi 
spectroscopy. Rabi oscillations cause changes of the 
qubit's magnetic moment and, therefore, excite the tank. 
The tank circuit accumulates photons which were emit- 
ted by the qubit. This approach is similar to the one 
in entanglement experiments with Rydberg atoms and 
microwave photons in a cavity. 

Indeed, quantitative information can be extracted from 
the noise spectrum S'y(ti') of the voltage fluctuations 
(the Fourier transform of the correlator Mv{t,t') = 

((l/2)[VT(i), VT(i')]+>) in the tank^ which incorporates 
the noise spectrum of the tank, Syr, supplemented by 
the qubit's contribution Svq,Sv = Svt{^) + SvQi(^), 



where 



(19) 



The qubit's part of voltage noise can be found from the 
stochastic equation l|ll() for the tank voltage: 



5'vq(cj) — X 



,2\2 



(20) 



Here Szzii-^) is a Fourier transform of the correlator 
{l/2[azfi{t),azfi{t')]_^), which describes internal fluctu- 
ations in the qubit (not related to the tank). Hand in 
hand with the tank's damping rate, 7t H15|I . the reso- 
nance frequency of the tank, Qt, is also shifted because 
of the qubit-tank interaction. 



= WT Jl - k^Lgl^x'zzi^r)- 



(21) 



The spectrum of voltage fluctuations has a peak near 
the resonant frequency of the tank lot, and, therefore, 
it contains information about a low-frequency compo- 
nent Szzi^^r) of the qubit spectrum. The equilibrium 
p art of th is spectrum peaks at the energy splitting uJc — 
vA^-Te^ of the tunneling doublet which differs signif- 
icantly from the frequency of the tank, ujc 3> lot- Be- 
cause of this, a contribution of equilibrium qubit noise 
to the voltage spectrum of the tank is expected to be 
negligibly small. An external microwave source with a 
frequency loq tuned in resonance with the energy split- 
ting of the qubit ujc, induces periodical variations of a 
population difference between the excited and ground 
state of t he qubit, which ar e characterized by a frequency 
Hr = \/ {A/uicYF'^ + 5q depending on the amphtude 
F of the microwave source as well as on the detuning 
Sq = loq — LOc- With non-zero bias, e ^ 0, left and right 
wells of the qubit potential have different energies. As a 
consequence, Rabi oscillations between the energy eigen- 
states will be accompanied by low-frequency transitions 
of the qubit from the left to the right well and back. The 
tank detects this kind of low-frequency noise which is de- 
scribed by the Lorentzian spectrum centered at the Rabi 
frequency flji with a linewidth dependent on the qubit 
decoherence rate F. Both the tank (F^) and the inter- 
nal heat bath (Fq) contribute to the decoherence rate 
F, F = Fo + Ft- It should be emphasized that the ex- 
ternal microwave field affects the qubit-bath coupling24 
that results in the distinction of the nonequlibrium decay 
rate F from its equilibrium counterparts F^: , Fj, entering 
the Bloch equations H16|) . 

An informative part of the spectrum of voltage fluctu- 
ations, 5'yQ(cL'), incorporates the qubit Lorentzian mul- 
tiplied by the transmission function of the tank having a 
sharp peak at the frequency lot '■ 



C't 



(^2 - Cj2)2 + ^2^2 

]hi 

(w2 _ f]|)2 + ^2r2 ■ 



(22) 
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FIG. 5: The spectral noise amplitude of the tank voltage for 
UHF powers Pa < Pb < Pc at 868 MHz. The bottom curve 
corresponds to the background noise without an HP signal. 
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FIG. 6: (a) Comparing the data to the theoretical Lorentzian. 
The fitting parameter is g ~ 0.02. Letters in the picture cor- 
respond to those in Fig. |S] (b) The Rabi frequency extracted 
from (a) vs the applied UHF amplitude. The straight line 
is the predicted dependence uin/uJT = \/ P/Po- The good 
agreement provides strong evidence for Rabi oscillations. 



The linewidth of the tank is assumed to be much less 
than the qubit's damping rate, jt ^ T. Because of this, 
as a function of frequency lo, the spectrum of vohage 
noise (|22ll represents a Lorentzian with a width and 
an amplitude which is given by a Lorentzian function of 
the Rabi frequency having the maximum near lot and the 
width F. Measurements of the noise spectrum amplitude 
at different values of the microwave power P allow to ex- 
tract information not only about the existence of Rabi os- 
cillations, but also about the nonequilibrium decoherence 
rate T of the qubit. We note that due to strong nonequi- 
librium conditions the populations of the qubit's levels 
are practically equal, and the noise spectrum amplitude 
does not depend on the temperature. The signal-to-noise 
ratio, 

peaks when fi^j = lot- At the same point, the backaction 
of the measuring device (tank) on the quantum bit which 
is described by the damping rate Ft, 

reaches its maximum as well. However, the tank con- 
tribution to the qubit decoherence drastically decreases 
with small detuning of the Rabi frequency from 
lutJt ^ \^R — ^t\ < T- At the same time, the efficiency 
of measurements, {SvQiuj)/SvT{uj))\i^=^rp, remains prac- 
tically unchanged. Since F <^ ^r,, the spectroscopic 
monitoring of Rabi oscillations with the low-frequency 
tank circuit falls into the category of weak continuous 
quantum measurements. 

The measurement setup as well as sample fabrication 
were similar to the ones described in the previous section. 



Microwave irradiation (UHF signal) was introduced to 
the sample by a commercial coaxial cable (between room 
temperature and ~ 2 K) and the resistive coaxial cable, 
known as ThermoCoax (between ^ 2 K and 10 mK). In 
order to reduce external disturbances, a 20 dB commer- 
cial attenuator was installed at 2 K. To measure Sy, we 
tuned the UHF signal in resonance with the qubit level 
separation. We found noticeable output signal only when 
ti;jjp/27r = 868 ± 2MHz, in agreement with the estimated 
splitting A/h ~ IGHz. Note that there is a difference 
of two orders of magnitude between co^^ and the readout 
frequency lot- Together with the high Qt, this ensures 
that the signal can only be due to resonant transitions 
in the qubit itself. This was verified by measuring Sy 
when biasing the qubit away from degeneracy. A signal 
exceeding the background, that is, emission of ~ 6 MHz 
photons by the qubit in response to a resonant UHF field 
in agreement with Eq. H24|l . was only detected when the 
qubit states were almost degenerate (cf. below Eq. if^ l. 
The measurements were carried out at nominal temper- 
ature T = 10 mK. No effect of radiation was observed 
above 40 mK (with 40 mK/ft-fce ~ 830 MHz, i.e. close 
to A/h) We plotted S'v(w) for different HF powers P 
in Fig. |5l As P is increased, grows and passes lot, 
leading to a non-monotonic dependence of the maximum 
signal on P in agreement with the above picture. This 
and the sharp dependence on the tuning of lo^j^ to the 
qubit frequency confirm that the effect is due to Rabi 
oscillations. 

For a quantitative comparison between theory and 
experiment, we subtracted the measured signal with- 
out an HF power from the observed Sy, yielding the 
qubit's contribution Svq = Sy — SyT{(^)- Subsequently, 
we extract ed th e peak values vs. UHF amplitude, 
<S'yQ,max(\/-P/A)) = max^Syqiuj) w SyqiLOr), where 
Pq is the power causing the maximum response; see 
Fig-Et- In the same figure, we plot the theoretical curve 
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for SvQ, max normalized to its maximum 5*0, 



SvQ,miLx{w) 

So 



2 9 



(5/2)^ 



(25) 

w = firi/cjT (= \/P/Po theoretically) and g = F/wt- 
The best fit is found for T « 0.02 lut ^ 8-10^ s'^ 
Thus, the life-time of the Rabi oscillations is at least 
TRabi — 2/r w 2.5 ^s, leading to an effective quality fac- 
tor QRabi = A/(^r) ^ 7000. These values substantially 
exceed those obtained recently for a modified 3JJ qubit 
(TRabi ~ 150ns) which is not surprising. In our setup 
the qubit is read out not with a dissipative DC- SQUID, 
but with a high-quality resonant tank. The latter is 
weakly coupled to the qubit (/c^ ~ 10~'^), 
the noise leakage to itm^ 



suppressmg 



IV. RESONANT SPECTROSCOPY 

In this section we show that the IMT can be also used 
for resonant spectroscopy, which is a well-known exper- 
imental method for investigation of quantum systems. 
As an example of such IMT application let us consider 
an interferometer-type charge qubiti^SiSLSS The device's 
core element is a single-Cooper-pair transistor - a small 
island, separated by two mesoscopic Josephson junctions, 
which is capacitively coupled to the gate. The transistor 
can be described by the Hamiltonian matrisSSiSS 



{Sn 



(26) 

where N is number of Cooper pairs on the island, 6n.m 
is Kronecker symbol, Ec — e^/2Cs is the single-electron 
charging energy expressed through the total island capac- 
itance Ce. The dimensionless parameter Ug ^ CgVg/2e 
is continuously controllable by the gate voltage Vg via 
the capacitance Cg. The effective Josephson energy 



ej{^) = [Ej^ + E% + 2EjiEj2 cos^] 



1/2 



(27) 



is a function of the total phase difference across both 
junctions (p = ipi +4>2, where Eji j2 and (pi^2 are Joseph- 
son coupling energies and phase differences of the first 
and second junction, respectively. 

If the transistor is closed by a superconducting loop 
with low inductance Lg, the total phase difference is 
(f) w 27r$e/$o and the ground-state curvature d?E^/d^\ 
can be obtained finding smallest eigenvalue of the Hamil- 
tonian matrix (|26|l as a function of (f>e. Using Q, we can 
calculte the phase shift of the tank inductively coupled 
to the charge qubit and compare it with experimental re- 
sults obtained by IMTi^ The principle of resonant spec- 
troscopy is very simple. If the qubit is irradiated by mi- 
crowaves with frequency coresponding to the energy gap 
between ground (n=0) and upper level {tl=\) the latter 
one becomes also populated. In this sense the microwave 



-0.2 - 



T3 



-0.3 - 



-0.4 * 



-0.5 





FIG. 7: Tank phase shift x vs. gate parameter Ug without 
microwave power (lowest curve) and with microwave power 
at different excitation frequencies. The data correspond to 
uJuHF/2n = 8.9, 7.5, 6.0 GHz (from top to bottom). The 
magnetic flux $e = $o/2 threading the interferometer loop 
provides a total phase difference S = n across the single- 
Cooperpair transistor. (For clarity, the upper curves are 
shifted.) 




FIG. 8: Energy gap A between the ground and upper bands 
of the transistor determined from the experimental data for 
the case S — tt. Some examples of these data are shown in 
Fig. |7| Dots represent the experimental data, the solid line 
corresponds to the fit (cf. text). 



irradiation acts like temperature i.e. suppresses the tank 
phase shift (see Ea ll8|l . 

Similar to the phase qubits, the interferometer- type 
charge qubit was fabricated out of Al by the conven- 
tional shadow evaporation technique, and was placed in 
the middle of the Nb coil by making use of a flip-chip con- 
figuration. The geometric loop inductance of the inter- 
ferometer was calculated to be Lg = 0.8 nH. The layout 
size of the junctions was 140 nm x 180 nm. Deviations 
from the nominal dimensions caused by the fabrication 
process were estimated from the micrograph of the real 
structure and found to be less than 15 %. The charg- 
ing energy was overestimated within the framework of 
the plate capacitor model from the junctions delivering 



Ec — 7 GHz. In fact and also in accordance with the 
experimental results below, this value is reduced due to 
the strong tunneling regimei^ The measurements were 
performed at mixing chamber temperature of 10 mK. 

The presence of the microwave power significantly 
changes the obtained dependence, namely peaks appear 
in the x(fT.g) curve (see upper curves in Fig.Cj). The peak 
position depends on the microwave frequency and does 
not depend on the amplitude(shape slightly depends). 
These peaks disappear when the phase bias is far from tt 
as well as at higher temperatures. Therefore, we believe, 
that they correspond to the excitation of the system from 
the ground to the upper state. 

The microwave induced transition (both the frequency 
of the microwave and the phase difference across the tran- 
sistor if = t: are fixed) from the ground to the upper state 
occurs only at certain value of the gate charge. From the 
position of the peaks on the x('^s) curves at different 
frequencies of the microwave, we have reconstructed the 
energy difference between ground and upper states as a 
function of the quasicharge on the island. The obtained 
dependence is shown in Fig. |H1 We fitted the experimen- 
tal data by using the numerical solution of the energy 
spectrum of the Hamiltoniani matrix 12611 . The best fit- 
ting parameters were found to be sj{tt) — 4.4 GHz and 
Ec = 2.2 GHz. This value of the Josephson coupling en- 
ergy is in very good agreement with the estimated value, 
and, as expected, the charging energy is smaller than 
estimated. 



NONRESONANT SPECTROSCOPY OF TWO 
COUPLED QUBITS 
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FIG. 9: Micrograph of the two-qubit system coupled to a 
resonant tank circuit. 




^J^Q (0.005 4>g per division) 



After the successful demonstration of quantum coher- 
ence in many types of superconducting qubits an obser- 
vation of entangled states in two coupled qubits presents 
the next step on the road to the quantum processor. 
The entangled states were recently observed in both the 
charge^ and the current-biased Josephson junction'^'* 
qubits. In this section we demonstrate that entangled 
states in a system of two inductively coupled flux qubits^ 
can be detected by the IMT,^ 

The system of two flux Al qubits inductively coupled 
to each other and to the Nb tank is shown in Fig. |51 
The area of each qubit and self-inductance, and critical 
current were Sq = 80 fim^, Lq = 39 pH, Ic ~ 400 nA, 
and Ec ~ 3.2 GHz, respectively. The mutual inductance 
between the qubits Mab = 2.7 pH was estimated nu- 
merically from the electron micrograph. The magnetic 
flux through the qubits was created by the dc compo- 
nent of the current in the coil Idci and by the bias cur- 
rent Idc2 through a wire close to one of the qubits. This 
allowed independent control of the bias in each qubit. 
The system of Fig. is described by the Hamiltonian 
H =^ Hq + Ht + Hint + Hdiss , where the two-qubit Hamil- 
tonian in the two-state approximation is expressed as'* 



FIG. 10; Normalized tangent of the current-voltage angle x 
in the tank vs. external flux bias <l>e at 50 mK. A relative 
flux bias f shift between the qubits is created by changing 
the current Idc2 in the additional wire. The shifted curves 
correspond to Idc2 ~ 27.3 fiA, while the central curve is for 
Idc2 ~ -2.7 fiA. The experimental data are denoted by the 
dots {Idc2 = 27.3 /iA) and triangles {Idc2 = -2.7 A'A). Sohd 
curves correspond to theoretical fitting. 



Ht is the tank Hamiltonian (a harmonic oscillator), the 
qubit -tank interaction is 



(29) 



and Hdiss describes the standard weak coupling of the 
qubits to a dissipative bath^. 



Here the coefficients are A, 



a/b 



Ma/b,Tla/b, where 



Hn 



- A.ai'^' ~ Abai'l -l-e.ai'^' +ebaf ) -I- Jdi'^) af^ , (28) 



Ma/b,T is the qubit-tank mutual inductance, La/b is the 
self-inductance and la/b is the amplitude of the persis- 
tent current in the corresponding qubit. In the stan- 
dard two-state approximation, the qubit current opera- 
tors are la/b = lajb'^'f^^^^ ■ The qubit biases are given by 

(a /a$o(/a; - 0.5 -|- /shift), Cfc = h^Q^jx " 0.5 ■\- ?7/shift), 

where the dimensionless flux ~ idc\ describes the held 
created by the niobium coil in both qubits, while the pa- 
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rameters /shift ~ Idc2 and 77 = Mtw/Maw < 1 give the 
bias difference between the qubits created by the addi- 
tional wire. Here Maw (Mtw) are the mutual inductances 
between the a (b) qubit and the additional dc wire (for 
our sample, Maw and Al^w were calculated numerically, 
yielding rj = 0.32). The qubit-qubit coupling constant 
J = Mablah is positive because the two qubits are in the 
same plane side to side, leading to antiferromagnetic cou- 
pling (according to the north-to-south attraction law). 

The application of the IMT for spectroscopy of two 
coupled qubits, similar to single-qubit problem (see sec- 
tion II) , can be conveniently discussed in terms of their 
magnetic susceptibility Xzz- In the linear-response ap- 
proximation the magnetic susceptibility of the two-qubit 
system Xzzi^^) is expressed through retarded Green's 
functions of the qubit operators cri"^'''', calculated with 
the equilibrium density matrix p = e*-^"^"^/-^ with Hq as 
in Eq. 1)28(1 . It can be generally assumed that the latter's 
eigenvalues E^, /i = 1,2,3,4, are non-degenerate and 
the eigenstates are orthonormalized, (j^|^) = Tak- 
ing into account the qubits' interaction with a dissipative 
environment ^ii^ we derive 



tanx 



■E 

tj.<i> 



E,y — E,, 



P„ 



(30) 



(31) 



where — eKp(—E^/T)/Z is a thermal population of 
the /Lt energy level, Z = Y^^eyi\){—E^/T), F^^ are deco- 
herence rates of the double-qubit system, and 

+ A.Afc(Ai|a?V>(Hf^i"V)- (32) 

At low frequencies w = ojt ^ l^-/^ — Ejj\ and in a weak 
damping regime, F^^ ^ \E^ — Ei,\, the decoherence rates 
F^,y have no effect on tanx, but are responsible for an 
equilibrium distribution in the system. 

The first two terms in Eq. 1(32(1 arc non-zero even if 
the two-qubit states are factorized. The first (second) 
term corresponds to the contribution of qubit a (h) and 
is nonzero near the qubit 's degeneracy point. These 
contributions are practically independent of whether the 
qubits' degeneracy points coincide or not. 

The second and third lines in Eq. I(32() describe co- 
herent flipping of both qubits, which is only possible for 
non-factorizable (entangled) eigenstates \v). There- 
fore the difference between the coinciding IMT dip of the 



two qubits and the sum of two single-qubit IMT dips 
provides a measure of how coherent is the two-qubit dy- 
namics (that is, whether entangled eigenstates of the two- 
qubit Hamiltonian Eq. 1(28(1 are formed) . This is a neces- 
sary condition for the system to be in an entangled (pure 
or mixed) state. 

The measurement results are shown in Fig. 1101 Com- 
parison of the single-qubit dips with the coincident IMT 
dip shows clearly that the contribution to tan x from the 
entangled eigenstates is significant. Indeed, the ampli- 
tude of the central dip in Fig.jTUIat T 50 mK is 1.12, 
compared to the sum of the single-qubit dips equal to 
1.69. It means that the entangled terms (second and 
third lines in Eq. I(32() ') are responsible for the negative 
contribution « —0.57 to tan^. 

At 50 mK the temperature is comparable to the char- 
acteristic energies in the two-qubit system (at the two- 
qubit degeneracy point the gap between the ground state 
and top excited state is ^ 100 mK). Since the character- 
istic measurement time in our approach is dictated by 
the much smaller tank frequency, lot, the system will 
have time to equilibrate. Indeed, an excellent quanti- 
tative agreement between the experiment (Fig. 110(1 and 
the theory Eq. I(31() confirms that the system is described 
by the equilibrium density matrix with the Hamiltonian 
Eq. ((28(1 (all the parameters of which we determined from 
the experiment). In other words, our system is in an equi- 
librium mixture of entangled two-qubit states. 



VI. SUMMARY 

We have shown that interferometer-type superconduct- 
ing qubits can be characterized by making use of the 
impedance measurement technique. Moreover, weak con- 
tinuous quantum measurements can be performed with 
this method. 
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